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We find concrete evidence for tlie presence of crisis-induced 
and Pomeau-Manneville Type-I intermittencies in an axisym- 
metric PDE mean-field dynamo model. These findings are of 
potential importance for two different reasons. Firstly, as far 
as we are aware, this is the first time detailed evidence has 
been produced for the occurrence of these types of intermit- 
tency for such deterministic PDE models. And secondly, de- 
spite the rather idealised nature of these models, the concrete 
evidence for the occurrence of more than one type of intermit- 
tency in such models makes it in principle possible that dif- 
ferent types of intermittency may occur in different solar-type 
stars or even in the same star over different epochs. In this 
way a multiple intermittency framework may turn out to be of 
importance in understanding the mechanisms responsible for 
grand-minima type behaviour in the Sun and solar-type stars 
and in particular in the interpretation of the corresponding 
observational and proxy evidence. 



I. INTRODUCTION 

Intermittency has been observed in a variety of real 
settings as well in a vast number of numerical models. A 
great deal of effort has therefore gone into understand- 
ing these modes of behaviour in the context of deter- 
ministic dynamical systems theory. These studies have 
demonstrated the existence of a number of different types 
of intermittency (such as Pomeau-Manncvillc [Q , Crisis 
0, On-off 1^ intermittencies), each with their own as- 
sociated signatures and scalings. Many of these forms 
of intermittency have in turn been concretely shown to 
be present in experiments and numerical studies of dy- 
namical systems in a variety of settings (see |H-p[ and 
references therein). 

An important potential domain of applicability of such 
behaviour arises in understanding the mechanisms under- 
lying the intermediate time scale variability in the Sun 
||7| - the occurrence of the so called Maunder or grand 
minima - during which solar activity (as deduced from 
the sunspot numbers) was greatly diminished ||7|,g| . This 
behaviour is also confirmed by evidence coming from the 
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analysis of proxy data 0. There is also some evidence 
for similar types of variability in solar-type stars |10[ . 

The idea that some type of dynamical intermittency 
may under pin the grand minima type variability in 
the sunspot record (the intermittency hypothesis [ll|) 
goes back at least to the late 1970's [||-|lj. This idea 
has been the subject of intense study over the recent 
years and has involved the employment of various classes 
of dynamo models, including ordinary differential equa- 
tions (ODE) (e.g. |l3yi5f| as well as partial differential 
equations (PDE) models (e.g. |l|-|l|,||l). In addition 
to the phenomenological evidence for the presence of 
intermittent-type behaviours in dynamo models [nq-EOJ, 
concrete evidence has recently been found for the pres- 
ence of particular types of intermittency in both ODE 
dynamo models |£l|,^ as well as a recently discovered 
generalisation of on-off intermittency, referred to as in- 
out intermittency |£^ , in PDE models |gj] . 

Here we wish to report concrete evidence for the oc- 
currence of two other types of intermittency, namely the 
crisis-induced and Pomeau-Manneville Type-I intermit- 
tencies, in PDE mean-field dynamo models. The organ- 
isation of the paper is as follows. In Sec. II we briefly 
introduce the model studied here. Sec. Ill summarises 
our evidence demonstrating the presence of these types 
of intermittencies in this model and finally in Sec. IV we 
draw our conclusions. 



II. MODEL 

Ideally one would wish to employ the full 3-dimensional 
dynamo models with a minimum number of approxi- 
mations and simplifying assumptions. Despite a num- 
ber of important recent attempts pB-p^, the difficulty 
of dealing with small scale turbulence makes a detailed 
and extensive self consistent study of such fully turbulent 
regimes in stars still computationally impractical (see e.g. 

In view of this an alternative approach in studies of 
stellar dynamos has been to employ mean-field models 
|l|,|6|jl|Jl|,|^,||. We should mention that there is an 
ongoing debate regarding the nature and realistic value 
of such models pO]. Nevertheless, 3-D turbulence simu- 
lations do seem to produce magnetic fields whose global 
properties (such as field parity and time dependence) are 
similar to those expected from corresponding mean-field 
dynamo models p3| . In this way mean-field dynamo 
models seem to reproduce certain features of the more 
complicated models and allow the study of certain global 



properties of magnetic fields in the Sun and solar-type 
stars (see for example [p3| , p4| )- 

The standard mean-field dynamo equation is given by 



dB 

'dt 



V X (u X B + q;B - ?7tV X B) 



(1) 



where B and u are the mean magnetic field and mean 
velocity respectively and the turbulent magnetic diffu- 
sivity rjt and the coefficient a arise from the correlation 
of small scale turbulent velocities and magnetic fields (a 
effect) ^^. We consider the usual algebraic form of a- 
quenching namely 



Qfo cos 9 



(2) 



where ao = constant and 6 is the co-latitude. 

We solve Eq. ^ in an axisymmetric configuration and 
in the following, as is customary [^ , we shall discuss the 
behaviour of the solutions by monitoring the total mag- 
netic energy, E = — !— /B^dl/, where /ip the induction 
constant, and the integral is taken over the dynamo re- 
gion. We split E into two parts, E — Ea + Es, where Ea 
and Es are respectively the energies of the antisymmetric 
and symmetric parts of the field with respect to the equa- 
tor. The overall parity P is given by P = [Es~Ea]/E, so 
P = —1 denotes an antisymmetric (dipole-like) pure par- 
ity solution and P = -1-1 a symmetric (quadrupole-like) 
pure parity solution. 

For the numerical results reported in the following sec- 
tion, we used a modified version of the axisynimetric dy- 
namo code of Brandenburg et al. (1989) [Q employed 
recently in pq]. These models are constructed from a 
complete sphere of radius R by removing an inner con- 
centric sphere of radius rg and a conical section of semi- 
angle ^0 about the rotation axis, from both the north and 
south polar regions (see |Q for details of the model and 
the relevant parameters). To test the robustness of the 
code we verified that no qualitative changes were pro- 
duced by employing a finer grid and different temporal 
step length (we used a grid size of 41 x 81 mesh points 
and a step length of 10~'^R'^/r]t in the results presented in 
this paper). For the following results we use Cq = —10*, 
which give the magnitude of the differential rotation and 
6o = 45°. The magnitude of the a-effect is given by the 
dynamo parameter Ca- 
in the next section we show in turn concrete evi- 
dence for the occurrence of crisis-induced and Pomeau- 
Manneville Type-I intermittencies. 



III. RESULTS 

A. Crisis— induced Intermittency 

As far as their detailed underlying mechanisms and 
temporal signatures are concerned, crises come in three 



varieties [|j. Here we shall be concerned with only one 
of these types, referred to as "attractor merging crisis" , 
whereby as a system parameter is varied, two or more 
chaotic attractors merge to form a single attractor. There 
are both experimental and numerical evidence for this 
type of intermittency (see for example (2|U] and references 
therein). In particular, this type of behaviour has been 
discovered in a 6-dimensional truncation of mean-field 
dynamo models pit]. 

Fig. |l| shows the plots of the energy and parity for 
the above model as a function of time, calculated with 
ro = 0.2 and Ca — 25.202 which show a bimodal be- 
haviour, switching intermittently between two different 
chaotic states. 
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FIG. 1. Example of crisis induced intermittency in a shell 
dynamo with a cut, with ro — 0.2, Ca ~ 25.202, Co — —10 
and So =45°. 

To determine the nature of this behaviour more pre- 
cisely, we have plotted in Fig. |2| the return maps for the 
PDE models (1), showing the attractors before and after 
the merging. As can be seen the resulting merged attrac- 
tor is, as expected, larger than the superposition of the 
two pre-existing attractors. 

These results can be taken as indications for the pres- 
ence of crisis-induced intermittency in this model. To 
substantiate this further, we recall that another impor- 
tant signature of this type of intermittency is the way r, 
the average time between switches, scales with the sys- 
tem parameter, in this case, Ca- According to Grebogi et 
al. g] , for a large class of dynamical systems this relation 
takes the form 



\Ca~C*a\ 



(3) 



where the real constant 7 is the critical exponent charac- 
teristic of the system under consideration and C* is the 
critical value of Ca at which the two chaotic attractors 
merge. 
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FIG. 2. Return maps showing the attractors in the PDE 
model (1) before (top panels) and after the merging (bottom 
panel). Note that as expected the merged attractor is larger 
than the superposition of the two previous attractors. 

The model under study here is a PDE system which is 
formally infinite dimensional. Such PDE models are nu- 
merically costly to integrate over long enough intervals of 
time (sometimes in excess of 5000 time units) necessary 
in order to obtain the scaling of the type (g). Further- 
more, the demonstration of such scaling requires a precise 
determination of the critical value C* which is difficult 
since as one approaches this value r diverges and the 
integration time becomes prohibitive. Despite these dif- 
ficulties, we have succeeded to obtain strong evidence for 
the presence of such a scaling as depicted in Fig. ^ with 
the corresponding 7 = 1.08 ± 0.05. Grebogi et al. B] 
conjecture that there may be a general tendency for 7 to 
be larger for higher-dimensional attractors. We do have 
a value of 7 higher than the previous one found for a re- 
lated six dimensional ODE dynamo model p^ but much 
lower than the value range suggested by Grebogi et al. 
Therefore, the conjectured range may need modification 
for large high-dimensional systems. 
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FIG. 3. Scaling of the average times between switches r 
as a function of {Ca — C%) for crisis induced intermittency for 
the model (§). The slope is found to be 7 = 1.08 ± 0.05 

There is also evidence for an enlargement of the final 
attractor after merging, as shown by the larger ampli- 
tudes of variation in the parity, in the sense that the 
parity gets closer to — 1 after the merging, as depicted in 
Fig. y. This helped us to numerically arrive at a better 
estimate for the critical value C*. 

These indicators, taken together, amount to strong ev- 
idence for the presence of crisis-induced intermittency for 
this model. 



B. Pomeau— Manneville Type-I Intermittency 

This type of intermittency, which is brought about 
through a tangent bifurcation, results in the system 
switching back and forth between a "ghost" periodic or- 
bit and sudden bursts of chaotic behaviour |l| . There are 
both experimental and numerical evidence for this type 
of intermittency (see for example HCTl and references 
therein). In particular this type of behaviour has been 
discovered in a 12-dimensional truncation of mean-field 
dynamo model p^ . 

To demonstrate the presence of this type of intermit- 
tency in the above PDE dynamo model, we have plotted 
in Fig. the energy and parity as a function of time for 
the parameter values tq — 0.7 and Cq = 28.0, which 
clearly demonstrates switches between nearly periodic 
behaviour and sudden bursts. We note that interestingly 
the energy in this case shows strong modulation which 
could be of interest in accounting for the occurrence of 
grand type minima in sunspot activity. 
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Example of Type-I intermittency in a shell dy- 
a cut, with ro = 0.7, C„ = 28.0, Ca = -10^ and ^ype-I mtermittency. 



FIG. 5. Power spectra of the time series in Fig. W for 



Another signature of this type of intermittency is pro- 
vided by the specific characteristics of its corresponding 
ower spectrum. By employing finite dimensional maps 
it has been shown that the corresponding spectra 
have a broad-band feature whose shape obeys approxi- 
mately the inverse-power law 1// for f > fs, where fs is 
the saturation frequency. Below this frequency there is a 
flat plateau induced by noise that causes arbitrarily long 
laminar phases to become finite. 

As further evidence for this type of intermittency in the 
model (1), we have plotted in Fig. o the power spectrum 
at Ca = 28.0, obtained by averaging over 16 different ini- 
tial conditions corresponding to different initial parities. 
As can be seen, the power spectrum shows both the flat 
plateau and the 1// power law scaling. 

Taken together, these indicators amount to strong evi- 
dence for the presence of Pomeau-Manneville Type-I in- 
termittency for this model. 



IV. CONCLUSION 

We have obtained concrete evidence, in terms of phase 
space signatures, spectra and scalings to demonstrate the 
presence of crisis-induced and the Pomeau-Manneville 
Type-I intermittencies in axisymmetric mean-field PDE 
dynamo models. Despite the rather idealised nature of 
these models, this is of potential importance since it 
shows the occurrence of two more types of intermittency 
(in addition to in-out intermittency recently discovered 
p4|) in these models which may in turn be taken as an 
indication that more than one type of intermittency may 
occur in solar and stellar dynamos. This suggests that 
any observational programme for identifying the mecha- 
nisms underlying grand minima type variability needs to 
take into account the possibility that multiple intermit- 
tency mechanisms may be operative in different stars of 
the similar type, or even in the same star over different 
epochs. This would also be of importance in the inter- 
pretation of proxy data. In this way a more appropriate 
hypothesis regarding such variability would be that of 
multiple-intermittency hypothesis. 
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